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We consider two different proposals to generate a time series with the same non-Poisson distri- 
bution of waiting times, to which we refer to as renewal and modulation. We show that, in spite of 
the apparent statistical equivalence, the two time series generate different physical effects. Renewal 
generates aging and anomalous scaling, while modulation yields no aging and either ordinary or 
anomalous diffusion, according to the prescription used for its generation. We argue, in fact, that 
the physical realization of modulation involves critical events, responsible for scaling. In conclu- 
sion, modulation rather than ruling out the action of critical events, sets the challenge for their 
identification. 

PACS numbers: 05.40.-a,05.65.+b,05.40.Fb 



I. INTRODUCTION 

The new field of complexity is attracting the attention 
of an increasing number of researchers, and it is trig- 
gering vivacious debates about its true meaning 0, 0, y| . 
Here we adopt the simple minded definition of complexity 
science, as the field of investigation of multi-component 
systems characterized by non-Poisson statistics. On intu- 
itive ground, this means that we trace back the deviation 
from the canonical form of equilibrium and relaxation to 
the breakdown of the conditions on which Boltzmann's 
view is based: short-range interaction, no memory and 
no cooperation. Thus, the deviation from the canonical 
form, which implies total randomness, is a measure of the 
system complexity. 

However, this definition of complexity does not touch 
the delicate problem of the origin of the departure from 
Poisson statistics. Here we limit ourselves to considering 
two different proposals, which we shall refer to as renewal 
and modulation. We shall show that these two proposals, 
although different, might lead to identical statistical re- 
sults, as far as waiting time distribution and correlation 
function are concerned. Thus, to a first sight, one might 
be tempted to conclude that they are indistinguishable, 
leaving no motivation whatsoever to prefer the one to the 
other. We shall prove that it is not so, and that there 
exist physical properties, whose observation allows us to 
distinguish the two proposals and to assess which one is 
correct for the given complex process under considera- 
tion. We shall prove, with numerical arguments, that, 
whereas renewal produces aging, modulation does not. 
This is the first result of the paper, confirming the theo- 
retical predictions of an earlier work [4| . We shall discuss 
also both modulation and renewal as source of diffusion 
generating fluctuations. We shall argue that the mod- 
ulation scaling depends, in principle, on the special way 
modulation is realized. We shall discuss two different pro- 



cedures to generate modulation and we shall show that 
both produce critical events. The critical events are re- 
newal events, leading a wealth of secondary events, and 
determining their Poisson properties. The critical events 
determine the scaling, which turns out to be anomalous 
only if the critical events are of non-Poisson type. This 
sets a challenge for the identification of the critical events, 
in general. 

The outline of the paper is as follows. In Section [D] 
we define the artificial time series under study in this pa- 
per, characterized by non-Poisson time distribution. In 
Sections IIIII and llVl we explain how to generate the non- 
Poisson waiting time distributions according to renewal 
and modulation prescriptions, respectively. In Section Ivl 
we discuss an aging experiment that turns out to afford 
an efficient criterion to distinguish between modulation 
and renewal. Section lVll is devoted to studying the diffu- 
sion process generated by either renewal or modulation, 
with the surprising conclusion that the physical realiza- 
tion of modulation generates unexpected renewal prop- 
erties. With Section IVlII we reiterate the importance of 
critical events, which, although ostensible in the case of 
renewal theory and invisible in the case of modulation, 
are responsible for anomalous scaling in both cases. 



II. TIME SERIES TO STUDY 

In this paper we shall discuss numerically and theoreti- 
cally the statistical properties of a time series {n} gener- 
ated from within two different approaches to complexity, 
renewal and modulation, in such a way as to produce in 
both ways the same time distribution with the following 
form: 
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where /i > 1. This distribution is properly normalized, 
and the parameter T, making this normalization possi- 
ble, gives information on the lapse of time necessary to 
reach the time asymptotic condition when ip{r) becomes 
identical to an inverse power law. The choice of the form 
of Eq. is dictated by the simplicity criterion. This 
form has been known for many years [f| , see for instance 
Ref. [6j, and following Metzler and Nonnenmacher |5j 
and Metzler and Klafter [3 we shall be referring to it as 
Nutting law. This form is also obtained by means of en- 
tropy maximization from a non-extensive form of entropy 
Q and, for this reason, is referred to by an increasing 
number of researchers as Tsallis distribution. 

We shall use the sequence {n} to distribute events on 
the time axis. The first event occurs at t = n, the second 
at time t = t\ + T2, and so on. The time intervals be- 
tween two consecutive events are called laminar regions, 
the reason being that the dynamical model here under 
study to illustrate renewal theory is an idealization of 
the celebrated Manneville map , with the time interval 
between two events representing the fluid regular state. 
We shall keep calling these time intervals laminar regions 
even when working with modulation, even if, in this case, 
as we shall see, the durations of different laminar regions 
are subtly correlated. 

For the theoretical discussion of this paper, it is useful 
to turn the sequence of times {t{\ into a diffusion gen- 
crating fluctuation Each time laminar region is as- 
signed through a coin tossing procedure either the value 
W or the value —W (W > 0), thus creating a dichoto- 
mous fluctuation £ (t) , which is interpreted as a stochastic 
velocity. We then shall study the time evolution of the 
coordinate x(t), obeying the dynamic prescription 
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We shall try to establish to what an extent the scaling 
of the resulting diffusion process depends on the origin 
of the fluctuation (modulation or renewal). This 
is done having in mind the interesting case of blinking 
quantum dots 0]> with W denoting the "light on" and 
— W the "light off" state. It is well known in fact that 
the waiting time distributions of "light on" and "light 
off" states depart from the Poisson condition. 

The central issue of the diffusion process under study 
in this paper has to do with scaling, namely, the property: 
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which is expected to hold true in the time asymptotic 
limit. The departure from ordinary statistical mechanics 
is signaled by either 5 ^ 0.5 or F(y) departing from the 
Gaussian form, or by both properties. We shall find the 
surprising fact that even with modulation the departure 
from ordinary statistical physics is determined by crucial 
renewal events, which seem to be an unavoidable conse- 
quence of any practical way we might adopt to realize 
modulation. 



III. RENEWAL 

As done in earlier work [TlllI^IT^IT^IT5] |. we are refer- 
ring ourselves to a simple dynamic process, as a prototype 
of renewal model. To make this paper as self-contained 
as possible, let us review here this simple model. Let us 
consider a particle moving within the interval / = (0, 1] 
driven by the following equation of motion 



dt V 
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and 



z>l, 



a > 0. 
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Due to the positivity of a the particle moves from the left 
to the right, and any time it reaches the border y = 1 
is injected back to a randomly chosen position yo, fitting 
the condition 



< y < 1. 



(7) 



The distribution density of sojourn times, ip{r), is evalu- 
ated as follows. First of all, we solve Eq. to determine 
the time necessary for the particle to reach the border 
moving from a given initial condition yo- This time is 
given by 
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The probability for the particle to get the border in the 
infinitesimal interval [r, t + dr] is determined by 



iP{ T )dT = p (y )dy Q . 



(9) 



We make the assumption of uniform back injection, which 
yields po(yo) = 1- Thus, we obtain the form of Eq. 
with 



(z-1) 



and 



T = 



(10) 
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It is interesting to notice that using Eq. leads to the 
following expression, for the mean sojourn time, (r) 



(r) = 



T 



Thus, using Eq. i|ll|) we express fi as follows 

T 
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This formula makes it evident that considering the case 
when the invariant distribution p eq {y) exists, and (t) < 
oo, yields > 2, which is the condition considered in this 
paper. 

In practice, we create first the sequence {yo(*)}, by 
means of a succession of random drawings of numbers 
from within the interval /. Then, using the transforma- 
tion of Eq. we associate yo(i) with Tj, thereby cre- 
ating the sequence {Ti}. We use this sequence to model, 
according to the renewal theory, the physical processes 
under discussion in this paper. Here we limit ourselves to 
using the sequence of 'laminar regions' {n} to distribute 
events on the time axis. As earlier stated, the first event 
occurs at t = ti, the second at time t = t\ + t-i , and so 
on. 

This model is inspired to the renewal theory 0] . The 
renewal character of the model is made evident by Eq. 
(|5|) . In fact, the values of yo are randomly chosen from 
a uniform distribution, < yo < 1. Any drawing does 
not have memory of the earlier drawings. Consequently, 
a laminar region does not have any memory of the earlier 
laminar regions. In literature there are many examples of 
renewal models yielding the distribution of Eq. Q • Here 
we quote first the model illustrated by Zaslavsky [l7| . 
We shall refer to this model as the hierarchical islands 
trap. It has to do with the condition of weak Hamil- 
tonian chaos, when regular islands are surrounded by a 
chaotic sea. The surface of separation is characterized 
by the island-around-island phenomenon: any iteration 
of a zooming process, yielding an increasing magnifica- 
tion, makes the same structures appear at smaller and 
smaller space and time scales. In other words, a parti- 
cle that reaches this separation surface through diffusion 
in the chaotic sea undergoes a stochastic motion that 
can be adequately described by a master equation cor- 
responding to an unbounded chain of states, \w)i, with 
i = 1, 2, • • • , oo. The first state of this chain, i — 1, is 
a doorway state, establishing a connection between the 
chaotic sea and the power law generating chain. The 
particle in this state can either jump back to the chaotic 
sea or forward to the second state of the chain, i = 2. 
From the second state the particle can either jump back 
to the doorway state, i = 1, or forward to the third state, 
z = 3. As the particle moves forward, namely i increases, 
the jump rate becomes smaller and smaller. Given the 
fact that both motion forward, increasing values of i, and 
backward, decreasing values of i, are possible, the parti- 
cle, sooner or later will jump back to the chaotic sea. A 
renormalization group approach jl7l | yields a distribution 
of waiting times corresponding, in the time asymptotic 
limit, to an inverse power law, and thus to a theoretical 
prediction compatible with the choice of Eq. £[]). The 
renewal nature of the process is ensured by the fact that, 
once the particle jumped back to the chaotic sea, the 
memory of the long sojourn in the power law generating 
chain is lost. 

In the specific case of blinking quantum dots, a model 
of the same kind as that of Ref. [lTj has been proposed 



by the authors of Ref. |3] • This model, in turn, is for- 
mally equivalent to that proposed years ago by Bouchaud 
[llil23, to explain the dynamics of glassy systems. These 
two models are significant examples of a wider category 
of renewal models, including the hierarchical islands trap 
model illustrated by Zaslavsky ■ It is straightforward 
to adapt the picture of hierarchical islands trap to the 
physics of blinking quantum dots. In this case, we say 
that the electron makes a jump from a state where spon- 
taneous emission of light is possible, to a state that can 
be thought of as the doorway state of a region where 
light emission is quenched. From this state the electron 
can jump either back or forward, the forward direction 
corresponding to a deeper and deeper embedding within 
the "light off" region. The forward and backward jump- 
ing rates become smaller and smaller with an increased 
embedding. The electron can come back to the doorway 
state and from there to the "light on" region. We can 
adopt an analogous model for the "ligth off" state. In the 
real case of blinking quantum dots this is a non-Poisson 
distribution, different from the non-Poisson distribution 
of the times of sojourn in the "light off" state. In this 
paper, for simplicity, we make the assumption that the 
two distribution are identical. This generates a sequel of 
times of alternate sojourn in the "light on" and "light off" 
state, with no correlation whatsoever among themselves. 
We shall see that the modulation approach apparently 
yields the same result, but it generates a subtle form of 
correlation among different sojourn times. 

IV. MODULATION THEORY 

We define as modulation theory any approach to non- 
Poisson distribution based on the modulation of Poisson 
processes. For instance, a double-well potential under the 
influence of white noise yields the Poisson distribution of 
the time of sojourn in the two wells [2l|. In the case of a 
symmetric double-well potential we have 

ip(t) = Aexp(-At). (14) 

The parameter A is determined by the Arrhcnius formula 

A = fccxp (^r)- (15) 

In the case when either the barrier intensity Q [2l| or 
temperature T [22j are slowly modulated, the resulting 
waiting time distribution becomes a superposition of in- 
finitely many exponentials. At least since the important 
work of Shlesinger and Hughes |2^, and probably ear- 
lier, it is known that a superposition of infinitely many 
exponentially decaying functions can generate an inverse 
power law. This, by itself, does not qualify the theory as 
modulation. It depends on the criterion adopted to gen- 
erate the sequence {Ti}, mentioned in Section^ If values 
of Ti , with different labels i are selected from different ex- 
ponential distribution, the resulting process is no doubt 



4 



renewal. To make it become a form of modulation theory, 
we have to draw a large sequence of time values, with the 
index i moving from i\ to i = i\ + Nd, with Nd >> 1, 
from the same Poisson distribution ip(r) — Aexp(— At). 

In recent times, the term superstatistics has been 
coined |^ to denote an approach to non-Poisson statis- 
tics, of any form, not only the Nutting (Tsallis) form, as 
in the original work of Beck 25]. We note that Cohen 
points out explicitly [24| that the time scale to change 
from a Poisson distribution to another must be much 
larger than the time scale of each Poisson process. Thus, 
we can qualify superstatistics as a form of modulation. 
Therefore, from now on we shall indifferently refer to 
this approach to complexity either as modulation or su- 
perstatistics. 

In conclusion, according to the modulation theory we 
write the waiting time distribution ip(t) under the follow- 
ing form 



function of the fluctuation £, 



V>(t) = / dAII(A)Aexp(-Af), 



(16) 

where 11(A) is the T distribution of order /i — 1 given by 
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This formula was proposed by Beck |25j and used in a 
later work pr|. 

We make the assumption of being able to generate time 
series with no computer time and computer memory lim- 
itation. Of course, this is an ideal condition, and in prac- 
tice we shall have to deal with the numerical limits of the 
mathematical recipe that we adopt here to understand 
modulation. The reader might imagine that we have a 
box with infinitely many labeled balls. The label of any 
ball is a given number A. There are many balls with the 
same A, so as to fit the probability density of Eq. (|17fl . 
We randomly draw the balls from the box and after read- 
ing the label we place the ball back in the box. Of course, 
this procedure implies that we are working with discrete 
rather than continuous numbers. However, we make the 
assumption that it is possible to freely increase the ball 
number so as to come arbitrarily close to the continuous 
probability density of Eq. (fT7|) . 

After creating the sequence {Xj}, we create the se- 
quence {ti\ with the following protocol. For any number 
Xj , the reader must imagine that we have available a box 
with another set of infinitely many balls. Each ball is la- 
beled with a number r, and in this case the distribution 
density is given by i^{t) = A exp (—At). 

To realize modulation, we adopt two different prescrip- 
tions: 

Prescription N.l We create a sequence {t^} by 
making Nd drawing from this box. Notice that accord- 
ing to the arguments of Cohen [24[ for modulation to be 
identified with superstatistics, it is necessary to make Nd 
very large, virtually infinite. Notice that the correlation 
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for any fixed A is equal to exp (—At). We note also that 
the smaller A the larger is the time interval corresponding 
to it. Consequently, for a proper definition of the effect 
of modulation on $^(t), we have to use the statistical 
weight n(A)/A [26j, which yields 
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This means that the waiting time distribution ip(r), of 
eq. (|l(jfl is proportional to the second time derivative 
of Eq. 119|) , this being a known consequence of renewal 
theory Hfl. Thus, for fi > 2, modulation and renewal not 
only yield the same iP(t), but also the same correlation 
function. 

As a consequence, studying the second moment of the 
diffusing variable x, through the average over the dis- 
tribution p{x 1 1) of the diffusing process is not sufficient 
to discriminate between renewal and modulation with a 
fixed number of drawings. We shall see, however, that 
the aging experiment of Section [3 allows us to distin- 
guish between modulation and renewal, in spite of the 
fact that, quite surprisingly, the adoption of this modu- 
lation prescription generates in the time asymptotic limit 
the same anomalous scaling as the renewal process. This 
is so because this kind of prescription generates by its 
own non-Poisson crucial events. These crucial events are 
imbedded in a sea of Poisson events that make them es- 
sentially invisible: the non-Poisson aging is annihilated 
by an overwhelming quantity of Poisson events that, yet, 
are not robust enough as to make the asymptotic scal- 
ing become ordinary. This delicate aspect is discussed in 
detail in Section IVT1 

Prescription N.2. In this case we keep drawing the 
waiting time numbers from the same Poisson distribution 
for a fixed amount of time T^. In this specific case 



$ 5 (t) = J dAil(A) exp(— At). 



(20) 



In this case the correlation function is obtained from the 
first time derivative of ip(t). Thus, in this case when 
(i > 2 the correlation function is integrable. The scaling 
of the diffusion process is expected to be that of standard 
diffusion, namely S = 0.5, while F{y) of Eq. © is a 
Gaussian function. 



V. AGING EFFECTS IN RENEWAL AND 
MODULATION THEORIES 

Aging of renewal processes is a term indicating the 
non-stationary nature of the probability distribution in 
the complete phase space cc, £, y, where the vector y 
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stands for the variables responsible for the dynamics of 
£, and hence responsible for ip(t). When starting from 
an off-equilibrium condition, the probability distribution 
of waiting times is affected by the evolution of the phase 
space distribution. As a consequence, the distribution of 
waiting times for the first event changes with the time 
t a , which is defined as the time delay between prepa- 
ration and observation. Let us explain this concept in 
some detail. We assume that preparation sets the sys- 
tem in a brand new, or "young" condition. This means 
a Gibbs ensemble with all the trajectories located at the 
beginning of a laminar region. Let us define t = as 
the time at which the observation of this system starts. 
The first measured waiting time is denoted by n. The 
first waiting time, at variance with the observation of the 
successive waiting times, does not necessarily correspond 
to the total time duration of a laminar region. In fact, 
the first laminar region could have started at some un- 
known time t — —t a , so that the first real waiting time 
is, let us say, given by t\ + t a . It should be noted that 
this is the key property generating aging and memory in 
the non-Poisson case. In fact, in the Poisson case the 
truncation of the laminar regions stemming from an ob- 
servation process with delay respect to preparation has 
the effect of altering both short and long-time laminar re- 
gions in the same way. As a consequence, a renewal pro- 
cess with Poisson statistics does not show aging. In the 
non-Poisson case, on the contrary, the short-time lami- 
nar region suffer from this truncation process more than 
the long-time laminar regions. Thus the distribution in- 
tensity at short times is lowered and at large times is 
enhanced, so generating a slower decay. The experimen- 
tal observation allows us to establish at which time the 
system was prepared, a form of long-standing memory 
that yet is compatible with the fact that the laminar re- 
gions showing up in the future do not have any memory 
of those that occurred in the past. In the model of Sec- 
tion IIIII this process of memory erasure is caused by the 
random back injection. 

In this section, we illustrate, with the help of a numer- 
ical example, how aging effects are deeply related to the 
renewal character of the process, whereas aging is annihi- 
lated by the slow modulation. Monte Carlo simulations of 
a renewal process with inverse-power-law distribution of 
waiting times, Eq. 1|T|). are performed using the approach 
proposed in Section fill! i.e., we generate a sequence {rj} 
by using Eq. JSJ. 

As far as modulation theory is concerned, the numeri- 
cal treatment is a little bit more delicate. The first step 
in the Monte Carlo simulation consists in generating a 
sequence of random numbers A in agreement with the T 
distribution 11(A) given in Eq. (|17|l . An approach sim- 
ilar to the one used in Section IIIII cannot be applied, 
because there is no simple analytical expression for the 
cumulative function related to T distributions. As a con- 
sequence, it is not possible to obtain a simple explicit 
expression for the function relating the random number 
A to a uniform random number in the interval / = [0, 1]. 



For this reason, we decided to use the rejection method 
suggested by Ref. j^. This method is based on the use 
of a majorant function /(A), i.e., a function that, for each 
point, takes values slightly greater than the correspond- 
ing values of the T distribution. The general idea is to 
draw random points in the plane, uniformly distributed 
under the graph of /(A). Then, the points under the 
graph of the T distribution are accepted and the others 
rejected. The random points are drawn as follows. First 
of all, we draw a random number A distributed as the 
probability distribution /(A) /A, where A > 1 is the area 
under the graph of /(A). This is easily done if /(A) is 
chosen with a simple analytical expression, in order to 
apply the method of Section IIIII Then, given A, a uni- 
form random number £ is drawn in the interval [0, /(A)]. 
The coordinates (A, Q define the random point. Finally, 
if £ is smaller than the corresponding value 11(A), then A 
is accepted, otherwise it is rejected. Without loss of gen- 
erality we fix T = 1 and we chose the following majorant 
function: 

(2i) 

According to Ref. , in order to obtain a certain num- 
ber of accepted numbers, rejection method requires a 
greater amount of random drawings, depending on the 
choice of the majorant function. We found that, given the 
majorant function described by Eq. I|21(l . with r] = 1.7 
and 7 = r/ ■ e (e is the Neper's number), the total number 
of drawings is about the double of the accepted draw- 
ings. Furthermore, we compared the histograms com- 
puted from the sequence of simulated random numbers 
with the relative probability distribution 11(A) given in 
Eq. i|17fl . with excellent agreement. Following Prescrip- 
tion N. 1, for each number A we have to draw Nd numbers 
from the exponential distribution described in Eq. I|14|) . 
representing a sequence of waiting times. This is eas- 
ily obtained with the same standard method described 
in Section IIIII Now, we used this approach to generate 
trajectories, i.e., artificial sequences of random waiting 
times. We compared the results for the aging analysis of 
trajectories characterized by the same exponent fi of the 
power-law, but generated from the renewal process and 
from modulation processes with different N^. 

We perform the aging analysis in the following way. 
Given the sequence of waiting times and an aging time 
t a , we compute the truncated waiting times, i.e., the dif- 
ference between each waiting time and t a . When the 
waiting time is shorter than t a , then we take the succes- 
sive waiting times until their sum exceeds t a . Then the 
truncated waiting time is defined as the difference be- 
tween this sum and t a . In this way, we obtain a sequence 
of truncated waiting times characterized by some aged 
probability distribution ip ta (t) . Aging is revealed if this 
distribution changes with t a . Note that the modulation 
process resting on Prescription N.l, for Nd — 1 is ex- 
pected to coincide with the prediction of renewal theory. 
To confirm this important property by means of numeri- 
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FIG. 1: Comparison between the function ^t tt of a re- 
newal process (continuous lines) and the of a modu- 
lation process with Nd = 1 (dots) at different values of 
t a — 0, 50, 100, 150, 200 (from the lower to the upper curves). 



cal calculation, in Fig. ^ we make a comparison between 
the renewal and the modulation process, with Nd = 1. 
Rather than measuring the function ip ta , we evaluate the 
aged survival probability ^t a , which is related to ipt a by 
the following relation: 

x/j ta ( T ')dr' = 1- J o 4h a {r')dr'. (22) 

As expected, the aged survival probabilities with the 
same t a , relative to the renewal and the modulation pro- 
cess respectively, coincide. In particular this happens for 
the brand new survival probabilities *&(t) = *&t a =o- The 
Survival Probability Functions ^t a ( T ) & H decrease with 
t and do not cross. Furthermore, for each fixed t, ^f ta (r) 
are increasing with t a . The comparison of Figs. [2 ED and 
01 reveals the lack of aging in the modulation theory with 
respect to renewal theory, which yields the maximum ag- 
ing effect for the smallest value Nd = 1, coinciding with 
the renewal condition. In the limiting case Nd — > oo the 
aging effect is annihilated. 



VI. DIFFUSION 

We address the problem of diffusion generated by mod- 
ulation using two distinct procedures: prescriptions N.l 
and N. 2 of Section lrTTl We shall use the continuous-time 
random walk to derive the asymptotic scaling of the pro- 
cess. The problem of diffusion with renewal has been the 
object of earlier work [2^| to which we refer the interested 
reader to appreciate the difference between renewal and 
modulation on this specific issue. Here we limit ourselves 
to noticing that the Levy scaling, namely 6 = — 1) 
and F of Eq. being a symmetric Levy (/i — l)-stable 
function, is a property of the central part of the spreading 
distribution, which is limited by ballistic peaks forcing it 
to become multi-scaling. 




FIG. 2: Comparison between the function ^ ta of a renewal 
process (dashed lines) and the function ^t a of a modula- 
tion process with Nd = 10 (dots), at different values of 
t a = 50, 100, 150, 200 (from the lower to the upper curves). 
The lowest curve represents the brand new function ^(t). 
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FIG. 3: The same as Fig. but with N d = 100. 
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FIG. 4: The same as Fig. H but with N d = 500. 



A. Continuous Time Random Walk: fixed number 
of drawings 

We now address the issue of revealing the crucial events 
corresponding to the use of prescription N. 1. At the 
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same time we plan to give a more accurate definition of 
crucial event, a property that so far has not yet clearly 
defined. The procedure that we plan to apply is as fol- 
lows. We study the diffusion process from t — to a 
given time t > 0, by assuming that for the particle to 
reach the position x at time t it is necessary to have n 
crucial events. By a crucial event we mean the drawing 
of a given value of A from the random distribution 11(A). 
Let us explain why the time of this drawing signals the 
occurrence of a crucial event. The concept of event im- 
plies unpredictability, a property shared by the random 
drawing of a given A and by the successive random 
drawings of the times r<'s from 



ip\(r) = Aexp(-Ar) 



(23) 



However, the property of being crucial is assigned only 
to the drawing of a given A. In fact, all the successive 
Nd time drawings adopt the same Poisson prescription, 
dictated by this leading drawing. It is remarkable that 
the sequence {n} would pass a random test, based on the 
ordinary correlation function, due to the random draw- 
ing. However, as we have seen in Section V, this sequence 
does not pass the aging test, which is a proof of the sub- 
tle correlation created by the persistent use of the same 
Poisson distribution for a large time span. 

We define as P(x, t) the probability of moving by a 
quantity x, either positive or negative, during time t, as 
an effect of the Nd drawings from the same Poisson dis- 
tribution. Analogously, with the symbol pW(i,i) we 
denote the probability that the particle moves by the 
quantity x, as a result of n crucial events, occurring in 
such a way as to ensure that the particle moves by the 
quantity x exactly in a time t. Note that after the occur- 
rence of these n crucial events, the particle is in general 
located at a distance x' from the origin different from 
x. With ri+l events the particle might overshoot that 
position. Thus the particle might make the remaining 
trip using the value of A drawn at the last crucial event, 
occurred at a time t' < t. In other words, we calculate 
p(x, t), i.e., the probability density of finding the particle 
in position x at time t using the formula 

p(x,t) = 

°° poo p-\-oo 

= J2 dt' p( n \x',t')N(x-x',t-t') 

„_1 JO J -co 



+N(x,t). 



(24) 



Here the quantity N(x, t) denotes the probability of mov- 
ing by a quantity x in a time t with no crucial event occur- 
ring. The Fourier-Laplace transform of p(x, t), p(k,u), is 
given by the expression 



p(k,u) 



1 



1 - P(k,u 



■N(k,u). 



(25) 



This expression is of crucial importance to determine the 
resulting scaling. To reach this important result we have 



to define first the function ip^(x,i). This is the proba- 
bility of moving by the quantity x in time t, with i time 
drawings from the same Poisson distribution. This im- 
portant function reads 



^{x,t) = / dAvi l) 0M)n(A), 



(26) 



where ip^\x,t) denotes the probability of moving by x 
in time i as a result of drawing i numbers from a given 
Poisson distribution, thereby corresponding to the same 
value of the parameter A. The functions tp^ (x,t) obey 
the following recursion relation: 

rt r+oo 

^\x,t)= dt' dx'i) { j- 1) (x',t')'il>X>(x-x',t-t / ), 

Jo J-oo 

(27) 



with 



Xexp(-Xt)- [S(x - Wt) + 5{x + Wt)} 



( 28 ) 

In conclusion, the physical meaning of the function 
ip^\x,t) of l|2t)|) corresponds to adopting the following 
procedure. With probability n(A) we draw a given A. 
Then, we draw i values of time intervals, and we move 
the particle by a quantity x with these i drawings. 

With these prescriptions we express the quantity 
P(x,t) = P<tt(x,t) as follows 



P (1) {x,t) = ^ Nd \x,t). 



(29) 



In other words, the transition by the quantity x in time 
t generated by a single crucial event is derived from Eq. 
(|27H by setting i = Nd- Let us remind the reader that 
our approach to modulation is based on drawing Nd times 
from the same Poisson distribution. 

Finally, we must express N(x, t). As earlier mentioned, 
the function N(x, t) denotes the probability of moving 
the particle by a quantity x in time t with no crucial 
event involved. Thus, 



N{x,t) 



dXn(X)N x (x,t), 



and 



N x (x,t) 



(30) 



(31) 



Nd nt f + OO 

V" / dt' I -/.-'.•" u t /' 

i=l 



dx'^- L> (x',t')y x (x-x',t-t'), 



where with ip^(x,t) we denote the initial condition 
5(x)S(t) In 13111 there are no crucial events involved be- 
cause the maximum numbers of time drawings is Nd — 1, 
not enough to generate a new crucial event, if we start 
at time t = with a crucial event. As a result of this 
partial number of drawings the particle reaches position 
x' in time t' . The remainder portion of space x — x' in 
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the remainder portion of time t — t' is traveled with no 
further time drawing, according to the prescription 



[S(x + Wt) + 5(x - Wt)] 



dt'e 



-At' 



(32) 

To establish the scaling produced by the modulation 
approach, we have to evaluate the Fourier-Laplace trans- 
form of ip^(x, t), i$!\(x,t) and P(x,t). Using the convo- 
lution theorem we have: 



i) {i) {k,u) 



+ 00 



dAn(A) 



X(u + X) 



and 



_{u + A) 2 + k 2 W 2 
(u + A) 



(u + A) 2 + k 2 W 2 ' 
From Eq. (|33|) we derive immediately 



P(k,u) 



+ OG 



dAn(A) 



X(u + A) 



(u + A) 2 + k 2 W 2 



(33) 



(34) 



(35) 



The evaluation of the asymptotic properties, u — » 
and k — > 0, is essential to assess the system scaling. If we 
make the assumption of sending u and k to first, and 
then integrating over A, we get 



p(k,u) = 



1 



u + k 2 W 2 ({) > 



(36) 



where the angle brackets with the subscript A indicate 
an average over the weight n(A)/A as in Eq. (|19Jl . Eq. 
yielding 



S = 0.5, 



(37) 



is equivalent to considering an ensemble of Poisson pro- 
cesses, each of them characterized by a fixed value of A, 
selected from the distribution 11(A) / A and then kept fixed 
forever. This suggests that the resulting scaling obeys the 
ordinary prescription S = 0.5. This is a correct property 
if /j. > 3. We note, however, that when /1 < 3, with \i > 2, 
due to the constraint of Eq. |Q , 



oc 

A ^0 



dX 



n(A) 

A 2 



(38) 



casting some doubts on the conclusion that S = 0.5. 

This is an incorrect conclusion indeed. In fact, with 
some algebra we find that 



X 2 texp(~Xt) 



x <^ [S(x - Wt) + 6(x + Wt)] + 



4 

9(Wt 



(39) 



N) 



Wt 

where denotes the step Heaviside function, and 

A 3 i 2 exp(-Ai) 



V4 3) (x,i) 



x <^ [S(x - Wt) + S(x + Wt)} + 3 



16 

9{Wt 



(40) 



N) 



wt 



The reader should compare these higher order distribu- 
tions with the expression of Eq. (|28|l . We see that in- 
creasing the number of drawing has the effect of weaken- 
ing the intensity of the ballistic component that, in the 
case of only one drawing, is known to generate anoma- 
lous diffusion [30l |. In fact, for only one drawing, using 
Eq. pgjl. we have 



(t) = * [S(x - Wt) + 5(x + Wt)] ^(t), (41) 



1 



where 



^(t) = (p-l) 



(t + Ty 



(42) 



Thus, with only one drawing, we obtain the same re- 
sult as that provided by the renewal theory. Notice that, 
with only one drawing, the dynamic process stems from 
the motion of two ballistic peaks (the two delta of Dirac 
of Eq. (28)). We see that already with Nd = 2, a contri- 
bution embedded between the to ballistic terms appear. 
We make a first estimation of scaling, based on neglecting 
these non-ballistic contributions. It is therefore evident 
that this prediction is not exact. 

To predict the effect of increasing Nd we notice that 
for t — > 00 the function ip^\x, t) will tend to a Brownian 
motion with a diffusion coefficient oc 1/A. This process, 
namely a renewal process with a single value of A is known 
to be described by the celebrated telegrapher's equation, 
whose solution is due to Cattaneo [3l| . If we define (t) 
the intensity of the ballistic peaks of tl)^(x,t), namely 
within a single value of A, and A^(£) the intensity of the 
peaks for the averaged distribution ^ l (a;,i), we find for 
the former the asymptotic expression 



A^(*./.) = i;<5 



XH l 



1! 



exp(— At). 
(43) 

In this expression we identify three factors with a spe- 
cific meaning. The factor 2~ l is a signature of the fact 
that for each ballistic contribution we tossed the coin i 
times, always getting the same result. The second term 
[6(x — Wt) + S(x + Wt)] means that we are consider- 
ing only the ballistic motion. Finally in the third term, 
A l i l_1 exp(— At)/(i — 1)! one can easily recognize the Pois- 
son distribution of obtaining the i-th success at the t-th 
trial, in a urn model. This indicates that the elapsed 
time is generated by i distinct drawings from the same 
Poisson distribution, with the same A. Let us average 
A! (je, t) over the weight 11(A) and write 



AW(i,t) 



1 



(i - 1)!2< 



[8{x - Wt) + S(x + Wt)] 



dAn(A)Aexp(Arf) 



(44) 
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From this form we obtain 

A w (x,t) = [S(x -Wt) + S(x + Wt)} (45) 

x (i-iy.&it + T)^- 1 ' 

The next step rests on replacing i with N d , thereby defin- 
ing the ballistic contribution to P(x, t), which is, as ear- 
lier pointed out, the probability of travelling by the quan- 
tity x, in the positive, x > 0, or negative, x < 0, di- 
rection, in time t. At this stage, we are ready to find 
the waiting time distribution of the times of sojourn be- 
tween two consecutive critical events, denoted by as as 
^ critical (t)- We obtain this waiting time distribution, 
concerning crucial events, by integrating A^ Nd ^ (x, t) over 
x, a procedure that yields 



critical 



(*) = 



l 



(N d -l)\2 N . 

jVt-ljJVi-l 



(t + Ty+ N i- 



— T (/x-l)M- ••(// + iV d - 2) 

(46) 



We see that the distribution of times of sojourn between 
two critical events reaches a maximum at 



{N d - 1)T 



(47) 



Thus, increasing N d has the effect of postponing the 
transition to the asymptotic regime. However, when the 
asymptotic regime is eventually reached, we see the emer- 
gence of an inverse power law with index /i. Due to the 
renewal character of the critical events, we get as a scal- 
ing S, the following expression 



S = 



1 



1' 



(48) 



which corresponds to a Levy diffusion process |29|. 

As earlier remarked, this prediction is not exact, given 
the fact that we have neglected the non-ballistic contri- 
butions to diffusion. We have considered only the delta 
of Dirac contributions, whose weight moreover tends to 
decrease upon increasing Nd- We wonder which might be 
the physical effect of the non-ballistic contributions that 
we have neglected. 

We estimate the role of the neglected terms by means 
of a heuristic approach, based on approximating the tele- 
grapher's diffusion with the ordinary Brownian diffusion. 
We denote this Brownian process with p\(x, t), express- 
ing the probability of traveling by the quantity a; in a time 
t, under the condition that a given A has been drawn at 
the beginning of the laminar region, at time t = 0, and no 
further drawing of A has occurred. We mentioned already 
that the diffusion coefficient is inversely proportional to 
A. ft is easy to find that 



Px(x,t) 



1 



{AttWH/X) 1 ' 2 



■ exp 



• A 



(49) 



To find the distribution M(x) for the space traveled in 
a time t c between two events, we assume t = t c ks N d /X 
and we average over A, namely, we write 



M(x) 



rpn- 



r( M 

A 



o {AirW 2 ^) 1 / 2 GXP \W 2 N d 



(- 



(50) 

X"- 2 exp (-AT) dX. 



Since the first moment (T cr iti C ai) of ip critical is finite, fol- 
lowing p9l | we assume that asymptotically, namely after 
a very large number of critical events, the central part 
of both P(x,t) and p(x,t) is approximated by a jump- 
ing process, with jumps of intensity M(x), occurring at 
regular times, namely the time interval between two con- 
secutive jumps, is fixed to be equal to {T cr iti C ai)- 

It is well known that if the second moment of M(x) is 
finite, the system falls in the basin of attraction of Gauss 
diffusion and 5 = 0.5. If, on the contrary, the second 
moment of x is divergent, the system falls in the basin of 
attraction of Levy diffusion, the diffusion process being 
alpha stable. This means that to recover finite moments 
we have to study the fractional moments with 
a < 2. Note that the index /i of the Levy inverse power 
law is fi = a+1. Thus, a heuristic way to establish if the 
process is of Levy kind or not, rests on the evaluation of 
the fractional moment (|a;| Q ), with the weight M(x) of 
Eq. I|5U|) . Then, for any a we must establish which is the 
threshold value of fi. By threshold value of fi we mean 
the value of n, below which the fractional moment (|a;| Q ) 
diverges. 

In other words, we set a < 2, and we look for the 
minimum value of ji ensuring 



dx\. 



dXX^ 



exp 



4W 2 N d 



- XT 



< oo. 



(51) 

We reverse the integration order and use the fact that 
the Gaussian moments are always finite, to get 



+oo 



dx\x\ a X a+i exp 



-(•?A) 2 
iW 2 N d 



= C. (52) 



Note that C is independent of of A. Thus the condition 
of EH becomes 



dXX h 



' exp (-AT) < oo. 



(53) 



Note that the only source of divergence for this integral 
is the limiting condition A — > 0. In fact the exponential 
factor exp(— AT) ensures convergence for all values of A. 
To ensure the convergence of this integral we must focus 
on the short values of A. We see that the integral conver- 
gence is ensured by /j, > a + 1. According to the earlier 
illustrated criterion, this yields for the scaling 6 = 1/a, 
S = !/(/!— 1), thereby confirming the Levy scaling of Eq. 
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B. Exact evaluation of asymptotic properties 

The Levy scaling emerging from a sea of Poisson 
events, which annihilate, as we have seen in Section V, 
any sign of renewal aging, is a disconcerting property. 
The arguments based on the analysis of M(x) afford al- 
ready a compelling proof of this surprising fact. However, 
as a form of double check, we plan to reach here the same 
conclusion through the study of p(k, 0), which, according 
to the Lew theory [32| . should bep(k,0) rj l — const-k a . 
Using (25) we get this asymptotic property provided that 
we show that P(k, 0) ~ 1 — const ■ k a . Note that we have 
to prove that a = fx — 1. For this purpose we focus on 
Eq. |j35J, we set u = 0, and we obtain: 



p(M) = 



r(M - 1) J 



y 



y- 



i 



/V 



dy. (54) 



At this point, we use the binomial expansion and the 
useful property: 



V 



N 



V 



2N 



Etc {t)v 21 



= i - 



X 1 ' N )y^ 



(55) 



Let us plug this expression into the integral, and divide 
it in two parts: 



F( M - 1) 



(56) 



y^ 2 e 



- Tk yd % 



^pW (N 



)y 



2i 



u-2 -Tky 



dy 



The first term gives T(/j, — l)/k^ and therefore the 
unity term in the zero-th order of P(k, 0). Looking at the 
second term, note that the integrand does not diverge in 
zero, at the infinity either: the integral can be considered 
as a Laplace transform (in Tk) of a finite quantity. So, 
at this first order, we obtain a constant. 

In conclusion, for the integral l|54|) we obtain the ex- 
pression 1 — const ■ fc p_1 , which is the same first order ex- 
pansion as that stemming from the renewal process, thus 
yielding H48|) . We think that at this stage the validity of 
Eq. (|48|l is widely proved. In spite of the approximation 
made to get the Levy scaling, Eq. I|48l) . this property is 
correct and is expected to show up in the time asymptotic 
limit. 



C. Continuous time random walk: fixed time Td for 
the action of the same A 

We have focused so far our attention on the modulation 
prescription N. 1. Let us consider now the modulation 
prescription N. 2. In the case where we keep drawings 
the times r from the same A for a time T d , assumed to 
be very large, we can apply the following argument. The 



probability for the particle to travel by the quantity x is 
given by 



M(x)=p(x,T d ) 



+ 00 



dxn(\) Px (x,T d ) 



(57) 



where p{x, T d ) denotes the probability of travelling by 
the quantity x, in the positive, x > 0, or negative, x < 0, 
direction, for a time X^, throughout which the random 
times Tj has been always drawn from the same Poisson 
distribution of Eq. I|2.'{ll . with a fixed value of A. Note 
that also with prescription N. 2 the crucial events corre- 
spond to the time when a new value of A is drawn. As a 
consequence of the random choice of A, as well as of the 
random choice of Tj, the function M(x) must be inter- 
preted as a distribution of totally uncorrelated numbers. 
Therefore, as done earlier, we have to focus now our at- 
tention on the fractional moments of M(x). Notice that, 
with prescription N. 2, the jumps occur at regular times, 
separated by the fixed distance T d , a property that makes 
easier and even more rigorous our scaling evaluation. 
For this prescription we have 



1 



{^w 2 T d /\y/ 2 exp \W 2 T d 

Substituting l(4^|l into equation ijBTjl. we obtain 



/ -x 2 X 



(58) 



Jo°° A 1 / 2 exp - AT) X^T^dX 



p ( X > Td) (4irW 2 T d y/ 2 T(p-l) 
Let us consider for simplicity T = 1: 

p{x, T d ) 



(59) 



T{n- 1) {^w 2 T d y/ 2 



x / exp 
/o 



m 2 T d 



1 A 



A^ 3/2 dA. (60) 



Now, using integral 3.381(4) of Ref.|H, 

f°° 1 
/ x v ~ x e-^ x dx = — 1», 
Jo V v 

with $l{u} > and KjV} > 0, we obtain 

r(/i-l/2) 1 ( AW 2 T d 



M(x) 



T(/i - 1) (4 7 rVF 2 T d ) 1 /2 \ x z + 4W 2 T d 



(61) 



(62) 

We see that for /i > 2, the second moment of M is 
finite. Thus, in this case the standard central limit the- 
orem applies, and consequently modulation yields 



5 = 0.5 



rather than 



M-l' 



(63) 



(64) 
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the scaling predicted by the renewal condition. Notice 
that in this case the correlation function ^ (£) is given 
by 



T 



T + t 



(65) 



This means that the correlation function is integrable. 
In Section VI D we shall illustrate two examples of dif- 
fusion equations based on the correlation function <E>£ (t) 
rather than on the critical events, visible or invisible. The 
adoption of the correlation function (f) is an important 
step of the ordinary procedures of statistical mechanics, 
which ignore the concept of critical event, and the con- 
cept of event as well. The emphasis of the conventional 
approaches is based in fact on the density time evolution, 
and the stochastic trajectories, with abrupt jumps is only 
indirectly taken into account. We shall see that modu- 
lation agrees with these theories with no critical events, 
only when diffusion becomes ordinary. Even in this case, 
however, resting only on the experimental evaluation of 
ip(t) would be misleading, since it would create the im- 
pression that non-ordinary statistical conditions apply, 
while the asymptotic behavior of the system is ordinary. 



D. Equations for processes with no crucial events 

We have seen that the physical realization of modula- 
tion implies the occurrence of critical events, these crit- 
ical events being determined by the procedure we adopt 
to change the parameter A with time. As an ideal exam- 
ple of equations without critical events we have in mind 
the following two generalized diffusion equations 

9 p(x,t) = W 2 f dt'$ c (t')-^p(x,t-t') (66) 



dt 



and 



^p(x, t) = W 2 j dt'<S> 6 (t')-^p(x, t) (67) 

Both equations are generated in such a way as to fulfill 
the dynamical constraint 



dt 



(x 2 (t))=D(t), 



where 



D(t) = W 2 [ dt'$ 6 (t'). 
Jo 



(68) 



(69) 



The latter equation is proven to yield a Gaussian form 
for p(x,t), whose width is given by (x 2 (t)) of Eq. I]68p. 
thereby yielding the scaling 



(3 



Note that this scaling is determined by the correlation 
function ${(£), which has the form 



T 



t + T 



(71) 



(70) 



with f3 — fi — 2or/i — 1, according to whether prescription 
N. 1 or N. 2 is used. 

It is interesting to notice that also the former equa- 
tion, Eq. yields the scaling of Eq. lj7U|) . as proved 
by the detailed analysis of Ref. 34] . It is convenient to 
notice that these two equations have a different origin, 
sharing only the assumption that no critical event exists. 
Both equations refer to the simple diffusion process de- 
scribed by dx/dt = £(t). Eq. (|67|l is exact under the 
assumption that £(t) is a Gaussian stochastic variable. 
For its derivation, see, for instance, Ref. [3^ and ref- 
erences therein. Eq. (|66|l refers to the case where the 
variable £ is a dichotomous process, and it is derived by 
using a Liouville-like approach ignoring the existence of 
trajectories with sporadic but abrupt jumps |34j . It is in- 
teresting to notice that both equations produce the same 
anomalous scaling of Eq. (|70JI even when the correlation 
function is not integrable. We have seen that the pre- 
scription N. 1, generates anomalous scaling, of the form 
of Levy, thus departing from the scaling of these theories, 
with the same correlation function, but with no critical 
events. With prescription N. 2 the critical events fall in 
the Gauss basin of attraction, the correlation function 
is integrable, and modulation produces the same scaling 
as the theories without critical events, but in this case 
the diffusion process is ordinary. This suggest that the 
emergence of anomalous scaling, of Levy kind, is always 
a signature of critical events, either visible or invisible. 



VII. CONCLUDING REMARKS 

Let us discuss which are the main results emerging 
from this paper. First of all, we found that the wait- 
ing time distribution is not, by itself, an unambiguous 
indicator of complexity. This confirms the conclusion of 
earlier work [36j from within a new perspective made at- 
tractive by the modulation, or superstatistics, approach 
to complexity [24l |25| . It is remarkable that the time se- 
quence generated according to the modulation prescrip- 
tion would pass the randomness test, based on the use of 
the correlation function method. In fact, as we have seen, 
the correlation enforced by modulation refers to the per- 
sistent use of the same Poisson prescription to generate 
waiting times that are otherwise totally uncorrelated. We 
have seen that the adoption of the aging experiment of 
Section V bypasses the limits of the conventional meth- 
ods of analysis, since it establishes beyond any doubt 
the difference between the renewal and the modulation 
character of a sequence of times. This is a fact of some 
physical importance for the physics of blinking quantum 
dots [ljj. The work of Brockmann et al. [37j has al- 
ready established aging in the intermittent fluorescence 
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of these new materials, and the more recent work of Ref. 
[38| confirms that this aging is due to the renewal charac- 
ter of the process. Thus, we conclude, in agreement with 
Refs. 0, H3) that the dynamic process responsible for 
intermittent fluorescence must be built up on the basis 
of a renewal perspective. 

Of course, the fact that superstatistics is ruled out as 
a proper perspective for the physics of blinking quan- 
tum dots does not mean, in any way, that there are no 
interesting complex processes in nature that might lay 
their roots on this attractive perspective. This remark 
is related to the second important result of this paper. 
Modulation, behind superstatitics, corresponds to a slow 
physical process whose realization does not rest on a 
unique procedure. However, it is a plausible that any 
realization of modulation might generate its own critical 
events. These crucial events are embedded in a Poisson 
sea, exerting a camouflage action, which makes it ex- 
tremely difficult, to reveal their existence. A few years 
ago the authors of Ref. [3|j established a technique of 
analysis of time series, through scaling evaluation, which 
was proved by later work j3(| to be an efficient method to 
reveal the statistical properties of critical events that, in 
spite of being invisible, namely, not even being recorded, 
determine the properties of the visible events, the events 
that are recorded. This paper proves that the physical 
realization of modulation rests on another class of invis- 
ible events, which determine the stochastic rather than 
the deterministic properties of the visible events. This 
sets the challenge for the detection of this new kind of 
invisible events. In Section VI we have seen that these 
invisible events establish the anomalous scaling of the 
process. On the basis of the earlier work of Ref. [2£j 
we think that this anomalous scaling refers to the central 
part of the probability distribution, and that the result- 
ing anomalous diffusion is multiscaling. There are good 



reasons to believe that with suitable developments along 
the lines of the work of Refs. [2^.l39l |. it will be possible 
to address successfully this challenging problem. This is 
left as a subject of future work. 

Finally, we want to remark that this paper establishes 
another important result, concerning the deep connection 
between critical events and anomalous scaling. There 
are theories that apparently produce complexity with- 
out involving critical events. In Section IVI Dl we have 
illustrated two diffusion equations generating anomalous 
scaling without involving critical events. Superstatistics 
is not compatible with any of them, but in the trivial case 
of ordinary diffusion. The reason is made clear by the 
results of this paper: the practical realization of super- 
statistics implies the occurrence of critical events, which 
might be either Poissonian or non Poissonian. In the sec- 
ond case a correspondence is generated between scaling 
and waiting time distribution, but the anomalous scal- 
ing is determined by these non-Poissonian critical events, 
rather than by the waiting time distribution. In fact, as 
proved in Section VI, the scaling is not determined by 
i/j(t), the visible waiting time distribution: the anomalous 
scaling is determined by ipcriticai (t) , namely the distribu- 
tion of the time distances between an invisible critical 
event and the next. This makes the role of critical events 
even more important, and it makes the renewal perspec- 
tive emerge also from within superstatitics, a theoretical 
perspective that seems to be apparently alternative to 
renewal theory. 
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